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Abstract For a (non-symmetric) strong Markov process X, consider the Feynman-Kac semi¬ 
group 

where A is a continuous additive functional of X associated with some signed measure. Under 
the assumption that X admits a transition probability density that possesses upper and lower 
bounds of certain type, we show that the kernel corresponding to Tf possesses the density 
pf(x,y) with respect to the Lebesgue measure and construct upper and lower bounds for 
pf{x, y). Some examples are provided. 
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1 Introduction 

Let (Xt)t>o be a Markov process with the state space K". For a Borel measurable 
function V : R” —R, we can define the functional At of X by 

At := [ ViX,)ds, t>0. (1.1) 

Jo 
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Suppose that lim(_>.o sup^ E^|At| = 0. Then, by the Khasminski lemma there exist 
constants C,b > 0, such that 


supE“el"^‘l < 


( 1 . 2 ) 


X 


see, for example, [11, Lemma 3] or [12, Lemma 3.3.7]. Estimate (1.2) allows us to 
dehne the operator 


r//(a;) :=E"[e^‘/(Xt)], x £ M", f > 0, 


(1.3) 


where the function / is bounded and Borel measurable. The family of operators 
(T/^)t>o forms a semigroup, called the Feynman-Kac semigroup. 

Feynman-Kac semigroup is well studied in the case of a Brownian motion (see 
[23, 24, 12, 3]); in particular, in [3] more general functionals are treated. The case of 
a general Markov process is much more complicated; see, however, [12, Chap. 3.3.2] 
and [24]. The essential condition on the process, stated in the papers cited, is that 
the Markov process X is symmetric and possesses a transition probability density 
Pt{x,y). 

In this paper, we construct and investigate the Feynman-Kac semigroups for a 
wider class of Markov processes. First, we construct the Feynman-Kac semigroup 
for a (non-symmetric) Markov process, admitting a transition density. We also treat 
a more general class of functionals At, that is, in our setting the functional At is not 
necessarily of the form (1.1), but is constructed by means of some measure va, which 
is in the Kato class with respect to the transition probability density of X (cf. (2.3)). 
The approach used in [8] allows us to show the existence of the kernel pf{x, y) of 
the semigroup {T^)t>o and to give its representation. The method from [8] relies on 
the construction of the Markov bridge density, which in turn employs the regularity 
properties of the transition probability density of the initial process X rather than its 
symmetry. 

In such a way, this prepares the base for the main result of the paper, which is 
devoted to the investigation of the Feynman-Kac semigroup for the particular class 
of processes constructed in [18]. In [20, 19], we develop the approach that allows us 
to relate to a pseudo-differential operator of certain type a Markov process possess¬ 
ing a transition probability density pt{x, y) and construct for this density two-sided 
estimates. In particular, such estimates provide an easily checkable condition when 
a measure vj belongs to the Kato class with respect to pt{x,y). This allows us to 
describe the respective continuous additive functional At and to show (1.2). Starting 
with the class of processes investigated in [18], we construct (see Theorem 3) the 
upper and lower estimates for the Feynman-Kac density pf{x, y). In particular, we 
show that the structure of such estimates is “inherited” from the structure of the es¬ 
timates on pt{x, y). In some cases when the upper bound on pt{x, y) can be written 
in a rather compact way, we can describe explicitly the Kato class of measures. For 
example, this is the case if pt{x, y) is comparable for small t with the density of a 
symmetric stable process; see also [4, Cor. 12] for refined results. In Proposition 4 we 
show that if the initial transition probability density possesses an upper bound of a 
rather simple (polynomial) form, this form is inherited by the Feynman-Kac density 
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Up to the author’s knowledge, in general, the results on two-sided estimates of 
p^{x, y) are yet unavailable. For X being an a-stable-like process, the estimates of 
the kernel pf‘{x, y) are obtained in [ 22 ]; see also [ 10 ] and the references therein for 
more recent results in this direction, including two-sided estimates on pf{x,y) in 
the case when the functional A is not necessarily continuous. The approach used in 
[22, 10] to construct the Feynman-Kac semigroup is based on the Dirichlet form 
technique. See also [5] for yet another approach to investigate Feynman-Kac semi¬ 
groups. 

The paper is organized as follows. In Section 2, we give the basic notions and 
introduce the main results. Proofs are given in Sections 3 and 4. In Section 5, we 
illustrate our results with examples. 

Notation 

For functions /, p, by / x 5 we mean that there exist some constants ci, C 2 > 0 
such that cif{x) < g{x) < 02/(x) for all x € R". By X ■ y and ||x|| we denote, 
respectively, the scalar product and the norm in R", and §" denotes the unit sphere in 
R". By i3b(R") we denote the family of bounded Borel functions on R". By C^(R") 
we denote the space of fc-times differentiable functions, with derivatives vanishing at 
infinity. By Ci, c and C we denote arbitrary positive constants. The symbols *, □, and 
() denote, respectively, the convolutions 

if*9)ix,y):= fix-z)g{z-y)dz, 

JR'^ 

ifng)ix,y):= fix-z)giz-y)widz), 

JR-^ 

and ^ 

if09)tix,y)-.= / ft-six,z)gsiz,y)widz)ds, 

Jo JR" 

where tJ7 is a (signed) measure. 

2 Settings and the main results 

Let X be a Markov process with the state space R". We call X a Feller process if the 
corresponding operator 

Ttfix)-.= ¥AfiXt) ( 2 . 1 ) 

maps the space Coo(R”) of continuous functions vanishing at inhnity into itself. As¬ 
sume that X possesses a transition probability density pt{x,y) which satishes the 
following assumption. 

PI. For fixed X G R", the mapping y 1—> Psix, y) is continuous for all s G (0, t], 
and the mapping s 1 —>■ psix, y) is continuous for all x,y G R". 

Recall some notions on the Kato class of measures and related continuous additive 
functionals. 

We say that a functional pt of a Markov process Xt is a PF-functional (see [13, 
§ 6 . 11 ]) if pt is a positive continuous additive functional, almost surely homogeneous. 
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and such that sup^, < oo. By additivity we mean that ipt satisfies the following 
equality; 

(fit+s = V’t +‘Ps o 9t, (2.2) 

where 0t is the shift operator, that is, Xg o 9t = Xt+g- The function Vt{x) := 
is called the characteristic of pt and determines pt in the unique way; see [13, 
Thm. 6.3]. 

A positive Borel measure w is said to belong to the Kato class Sk with respect 
topt{x,y)\f ^ 

lim sup / / pg{x,y)vu{dy)ds = 0. (2.3) 

*“*■0 xeK" Jo jR" 

By [13, Thm. 6.6], the condition zu G Sk implies that the function 



Pg{x,y)w{dy)ds 


(2.4) 


for which the mapping x Xt{x) is measurable for all t > 0, is the characteristic of 
some W -functional pt. 

Let vj = 07+ — vj~ be a signed measure such that G Sk with respect to 
Pt{x,y). Then 


xt-= / Ps{x,y)w^{dy)ds 
Jo JR" 


(2.5) 


are the characteristics of some VL-functionals A^, respectively, that is, there exist 
such that x^ix) = WAf. Since for such functionals we have 


lim supE^A^ = 0, 


then estimate (1.2) holds true, and thus the Feynman-Kac semigroup {T^)t>o for 
A := At - A,- is correctly defined. 

To show that the semigroup (T/^)t>o can be written as 

Tt'f{x)=[ f{y)pt{x,y)dy, 
jR" 

and to find the representation of the density pt{x, y) in terms of the probability den¬ 
sity of the initial process, recall some notions on Markov bridge measures. 

Denote by (St)t>o the admissible filtration related to X. A Markov bridge Xt’^ 
of Xt is a Markov processes conditioned by Xq — x and Xt — y. In the proof of 
[8, Thm. 1], it is shown that under PI there exists the corresponding Markov bridge 
measure y on Tt_ for {t,x,y) such that pt(x,y) > 0. We denote by E^ y the 
expectation with respect to P^ j,. 

The next proposition is essentially contained in [8, Thm. 1], but we reformulate 
the result in the way convenient for our purposes. 

Proposition 1. Let X be a Feller process, admitting the transition probability density 
Pt{x, y), for which assumption PI holds. Let w = w'^ — w~ be a signed Borel mea¬ 
sure, G Sk, ond At = Af — Af, where A^ are continuous additive functionals 
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with characteristics (2.5), respectively. Then 

Tt'^fix)= f f{y)pf{x,y)dy for any f £ Bb{W^), 

J {y.pt{x,y)>t)} 


where 

Pt{x,y)=pt{x,y)'¥}^^ye^\ x,yeW^,t>0. (2.6) 

Remark 1. When 26 is a Brownian motion, the statement of Proposition 1 is known, 
see [23] and also [3]. The construction from [3, 23] can be extended to the case of a 
symmetric Markov process, see [24]. On the contrary, the construction presented in 
[ 8 ] relies on PI and does not require the symmetry of the initial process. 

Proposition 1 implicitly gives the representation of the function pf{x, y). How¬ 
ever, when one wants to get quantitative information about pf{x,y), like the up¬ 
per bound on pf{x,y), estimation of the expectation in ( 2 . 6 ) appears to 

be non-trivial. Instead, for some class of Feller processes, we can use another ap¬ 
proach, which enables us to get explicitly an upper estimate of pf{x, y). Namely, in 
[18] we formulated the assumptions under which one can construct a Feller process 
possessing the transition probability density pt{x,y) satisfying assumption PI and 
admitting upper and lower bounds of certain form. In order to make the presentation 
self-contained, we quote this result below. 

Let 

Lf{x) := a(x)-V/(x)-f / (/(x-f m)-/( x)-wV/(x)l{||„||<i})m(a;,u)/r(du), 

(2.7) 

where / € and p is a Levy measure, that is, a Borel measure such that 

/ (||u||^ A l)/r(du) < cx). 

JR" 

Assume that p, satisfies the following assumption. 

Al. There exists /3 > 1 such that 

sup {r£) < j3 inf q^{r£) for all r > 0 large enough, 
reS" reS" 

where 

q^iO-=f [(?'w)^Al]/i(du), q^{£)-= f {^-ufpidu). ( 2 . 8 ) 

7r" 7|u-{|<i 

Assume that the functions a{x) and m(x, u) in (2.7) satisfy the assumptions A2- 
A4 given below. 

A2. The functions m{x, u) and a{x) are measurable, and satisfy with some con¬ 
stants hi, 62 , > 0 , the inequalities 

hi < m{x,u) < 62 , |q( 2 ;)| < b^, x,u G R". 
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A3. There exist constants 7 € (0,1] and 64 > 0 such that 

\m{x,u) — m{y,u)\ + ||a(a;) — a(y)|| < 64(||a; — t/P A l), u, x, y € ffi.". 

(2.9) 

A4. In the case /3 > 2, we assume that a(x) = 0 and the kernel m(x, u)y{du) is 
symmetric with respect to u for all x G K". 

Denote by /low and /up the functions of the form 

/iow(a;) := ai(l - a 2 ||x||)_^, fup{x) := x S K", (2.10) 

where > 0, 1 < i < 4, are some constants. 

Finally, define q*{r) := sup^ggn {r£), r > 0. It was shown in [17] (see 
also [20]) that condition A1 implies that 

q*{r) > r > 1. 

Note also that the continuity of q^ in ^ implies the continuity of q* in r. Therefore, 
we can define its generalized inverse 

pt :=M{r : q*{r) = 1/t}, fe(0,1]. (2.11) 

Theorem 2 ([18]). Under assumptions A1-A4, the operator {L, (R") extends to 
the generator of a Feller process, admitting a transition probability density ptix, y). 
This density is continuous in (t, x, y) G (0, 00) x K." x R", and there exist constants 
Oi > 0, 1 < z < 4, and a family of sub-probability measures {Qt, f > 0} such that 

PUiow{{x-y)pt) <Ptix,y) < Pt{fnp{Pf)*Qt)ix-y), t G (0,1], x,y G R”, 

(2.12) 

where /low and /up are functions of the form (2.10) with constants Qi, and pt is 
defined in (2.11). 

The constructed process is a Levy type process. In the “constant coefficient case,” 
that is, where a(x) = const and m{x,u) = const, (2.7) is just the representation 
of the generator of a Levy process; in other words, a Levy type process is the pro¬ 
cess with “locally independent increments.” It is known (cf. the Courrege-Waldenfels 
theorem, see [16, Thm. 4.5.21]) that if the class C'/°(R") of infinitely differentiable 
compactly supported functions belongs to the domain D{A) of the generator A of a 
Feller process, then on this set (R") the operator A coincides with L -f “Gaussian 
component.” Thus, the class of processes satisfying the conditions of Theorem 2 is 
rather wide. 

Let us show that, under the conditions of Theorem 2, we have 
Ptix, y) > 0 for all t > 0, x,y G R". 

We find the minimal N such that the distance from x to y can be covered by N balls 
of the radius smaller than (2a2Pt/jv)~^ (where 02 > 0 is the constant appearing in 
/low in (2.12)), that is, the minimal N for which 


N 


a2Pt/N 


(2.13) 
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Observe that q*{r) < cir^, r > 1, implying C 2 t < pt for all t small enough. 
Hence, (2.13) holds with N > 02C2lk-y||) ^ Therefore, putting = x and = y, 
we get 


Pt{x,y) = 


y) = ■■■ {Y\pt/N{y^-l,yi)]dyl...dyN 

jR" JR" J 


N 


> 


— ^OPtjNj 


' B(yo,i2a2Pt/N) '■) B(!/N‘-l,(2a2Pt/Ar) 


Y[pt/N{y^-l,yi)dy^ 


Nn 


where in the last line we used that 


Pt/N{yi-i,yi) >‘^ ^aiPt/N Morally, € B{yi_i,{2a2Pt/N) ^)- 

Thus, the transition probability density pt{x, y) is strictly positive. 

Finally, for a signed Borel measure zu, dehne 

/i(r) := sup |t77|{y : \\x - y\\ < r}, (2.14) 

X 

where |i37| := +'cu~ is the total variation of w. Denote by h the Laplace transform 

of h. 

The following theorem is the main result of the paper. Let to G (0,1] be small 
enough. 

Theorem 3. Let X be the Feller process constructed in Theorem 2. Take a signed 
Borel measure w such that its volume function (2.14) satisfies 

p^+^h{ps)ds < Ct‘^, t € [0,1], (2.15) 

with some constants C, C > 0, where pt is given by (2.11). Then 

a) There exists a continuous functional At such that 

E^At= [ [ ps{x,y)w{dy)ds] 

Jo JR" 



b) The semigroup {T^)t>o is well defined, and its kernel possesses a density 
p^{x, y) with respect to the Lebesgue measure on ffi."; 

c) There exist constants > 0, 1 < * < 4, and a family of sub-probability 
measures {3^t, f > 0} such that for t G (0, fo] <^nd x,y G M", 

Ptf\oy,{{x-y)pt) <ptix,y)<Pt{fupipf)*‘Jit){y-x); (2.16) 


here /low ond /up are the function of the fonn (2.10) with some constants at, 
1 < i < 4. 
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Remark 2. In general, Ui in estimate (2.16) are some constants, that may not coincide 
with those in estimate (2.12). In order to simplify the notation, we assume that in 
Theorem 2, oi = 03 = 1, 02 = a, and 04 = b. 

Assumption (2.15) can be relaxed, provided that more information about the ini¬ 
tial transition probability density is available. Put 




1 

f^(l+ ||x||/fl/a)d+a’ 


t> 0, X € K". 


(2.17) 


Note that for d = n, this function is equivalent to the transition probability density of 
a symmetric a-stable process in R." (that is, the process whose characteristic function 
is ). Denote by 'Xn,a the class of Borel signed measures such that 


lim sup 

t — ^0 7* 


\w\{y : ||x - y\\ < s} 

^n+1—a 


ds = 0. 


(2.18) 


The following lemma shows that for d > n — a the Kato class of measures with 
respect to gt(a; — y) coincides with Xn,a- The proof uses the idea from [4], and will 
be given in Appendix A. 

Lemma 1. A finite Borel signed measure vo belongs to Sk with respect to Qt{x — y), 
given by (2.17) with d > n — a, if and only if\w\ £ 'Xn,a- 

Corollary 1. In particular, it follows from Lemma 1 that w £ Sk with respect to the 
transition probabiility density of a symmetric a-stable process if an only ifzu £ 'Xn,a- 

In the proposition below, we state the “compact” upper bound for pf{x, y). 

Proposition 4. Let X be a Feller process satisfying the conditions of Proposition 1, 
and in addition assume that the transition density pt(x,y) of X is such that for all 
t £ (0,1], x,y £ K", the inequality 


Pt{x,y) <cQt{x-y), f e (0,1], x,?/G K", (2.19) 

where the function Qt{x) is defined in (2.17) with d > n — a. Suppose that w £ %n,a- 
Then 

pt{x,y) <CQt{x-y), f G (0,1], x,y G R”. (2.20) 

Remark 3. a) For X being a symmetric a-stable-like process, such a result is known, 
see [22]. In particular, the upper bound (2.20) holds with n = d. In our case, X is 
from a wider class; in particular, we do not assume the symmetry of the initial process, 
and the method of constructing the Feynman-Kac semigroup is completely different. 

b) In view of Lemma 1, under the assumptions of this proposition, we can take 
XU £ Xn,a rather than zu £ Sk vvith respect to gt, which is more convenient for 
usage. 

In Section 5, we provide examples that illustrate Theorem 3 and Proposition 4. 
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2. 1 Discussion and overview 

1. On continuous additive functionals. Loosely speaking, there are two approaches 
for constructing continuous additive functionals. One approach, which we de¬ 
scribed previously, relies on the Dynkin theory of kk-functionals. Another ap¬ 
proach, based on the Dirichlet form technique, establishes the one-to-one cor¬ 
respondence between the class of positive continuous additive functionals and 
the class of smooth measures, see [14, Lemmas 5.1.7,5.1.8] or [15, Thm. 5.1.4] 
in the case when the process under consideration is symmetric; see also [21, 
Thm. 2.4] for the non-symmetric case. In this paper, we use Dynkin’s approach 
as more appropriate in our situation, in particular, we do not assume that the 
initial Markov process X is symmetric. Our standard reference in this paper 
is [13]. 

2. On the generator of {Tf^)t>o. Suppose that the Markov process X and the 
positive functional At are as in Proposition 1. In this case, the semigroup 
(Tt^) t>o is contractive, and thus there exists a sub-Markov process with transi¬ 
tion sub-probability density pf{x, y). Formally, we can describe the generator 
of {Tt^)t>o as 

L'^ = L-zu, (2.21) 

where L is the generator of the semigroup associated with X, and vu is the 
measure appearing in the characteristic of At (cf. (2.4)), see [13, Thms. 9.5,9.6] 
for the (equivalent) formulation. Nevertheless, in this framework the problem 
of dehning the domain D{L^) of still remains open. In the general case, 
that is, when A can attain negative values, in order to dehne the generator of 
(non-contractive) semigroup {Tf^)t>o, we can use the quadratic form approach, 
see [1, 2], and also [9]. 


3 Proof of Theorem 3 

3.1 Proof of statements a) and b) 

a) By the upper bound in (2.12) on pt{x, y) (see also Remark 2), (2.15) implies that 
vj € Spc'. 


sup / / Ps{x,y)\w\{dy)d.‘s 

cgR”Jo JR" 

/ Ps/up((y -X- z)ps)Qs{dz)\w\{dy)ds 

n poo 

/ p'l\w\{y ■. \\y - X - z\\ <v/ps)e~^'"dvQs{dz)ds 
X Jo 


a;GR”J 0 JR" 

< sup / 
a;GR" JO 


<b Pg'^h{bps)ds —>-0, f —>■ 0. 


Hence, applying [13, Thm. 6.6], we derive the existence of a continuous functional 
At with claimed characteristic. 


Statement b) is already contained in Proposition 1 . 
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3.2 Outline of the proof of c) 

For the proof of Theorem 3(c), we use the Duhamel principle. First, we show that the 
function y) satishes the integral equation 

Pt{x,y)=ptix,y)+ / pt-s{x,z)pf{z,y)m{dz)ds, (3.1) 

Jo 

provided that the integral on the right-hand side converges. We show that if the series 


Trt{x,y) := '^pf^{x,y) (3.2) 

k=l 

converges, then it satishes Eq. (3.1). We derive an upper estimate for the convolutions 
pf^{x, y), which guarantees the absolute convergence of the series and allows to hnd 
the upper estimate for 7rt(a;, y). 

Second, we show that on (0, tg] x K" x R" the solution (3.2) to (3.1) is unique in 
the class of non-negative functions {f{t, x,y) > 0, t € (0, to], x,y € R"} such that 


f{t,x,y)dy < C 


for all t e (0, tf\, X € R". 


(3.3) 


We use the standard method, based on the Gronwall-Bellman inequality. 

Finally, observe that the kernel pf{x, y) of belongs to the class of functions 
satisfying (3.3). Indeed, since for At we have (1.2), it follows that 

|T/'/(a;)| < ciE^el^‘l < C 2 , / G ^^(R”), x e R", f G (0,fo]- (3.4) 

Thus, j/) = nt{x,y) on (0,fo] x R" x R”. 

Before we prove that (3.2) is the solution to Eq. (3.1) on (0,fo] x R" x R", 
let us discuss a simple case when w is the Lebesgue measure on R". In this case 
h{r) = c„r”, and thus assumption (2.15) is satisfied; 

[ p'^~^'^h{ps)ds = Cnt. 

Jo 

Therefore, the procedure of estimation of convolutions reduces to those treated in [ 18, 
Lemmas 3.1, 3.2]. 

Rewrite the upper bound in (2.12) as 


Ptix,y)<Cit *Qt){,y - x), 

(3.5) 

where Ci > 0 is some constant. 


gP{x) := t^/^gtix), 

(3.6) 

and (cf. Remark 2) 


gtix) :=p-/,p(ptx)=pre“'"‘'"'- 

(3.7) 
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This modification is technical, but proves to be useful for estimating the convolutions 
pf^{x,y). Let us estimate pf^{x,y). Take now a sequence {0k)k>i such that 0 < 
9k+i < Ok, 0i = 1 , and put 

:= t^^’^gtiOkx), A; > 1. (3.8) 

Since pt is monotone decreasing, for 0 < s < |, we have pt-s < Pt/ 2 - Note that 
Pt Pt/2’ follows from condition A1 and the definition of pp, see [20] for the 
detailed proof. Then, for 0 < s < f/2. 


^ / gt-s{0k-ix - 0k-iy)gs{0k-iy)dy 

^^k/2g-n^ f g^_^(^Q^_^x-y)gs{y)dy 

hl‘lf\—n f T) T) — ng ^_2,\-\.\z\) —bpsil — 77^— )\z\ , 

< t ' Qp-l / P^-sPs^ 


Jr-^ 

= Dkgl"\x), (3.9) 


where Dk = c{9k-i—0k)~^, c = ci e~^^^^dz, and in the second line from below, 
we used the triangle inequality and monotonicity of pt- In the case t/2 < s < t, 
calculation is similar. 

By induction we can get 


\pf^ix,y)\ <Ckt^^ *QP)iy - x), k>2, (3.10) 


where 


and for k > 2 




n 


nk/2) 


i=2 


QPidw) := f j (l-r)('= ^/‘^Qffl,^{dw-u)Qf}{du)dr. 


Since {Qi^\ t > 0, k > 1 } is the sequence of sub-probability measures and g[^'^ (x) < 
^n^fc/ 2 , .^g obtain 

jpf^(x,y)l<Ckt^-^p^. 

Thus, to show the absolute convergence of the series pf^ix, y), we may check 

that < oo. However, the behaviour of Cfc as A ^ oo is rather complicated. 

To see this, take, for example, = 5 + ^ - Then 


Ck = d 


,fc-i 


L'=(l/ 2 ) 

T(m 


(2'=A!(A-1)!)' 
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and thus Ck explodes as k ^ oo. Therefore, this procedure of estimation of con¬ 
volutions is too rough, and needs to be modified. For this, we change the estimation 
procedure after some finite number of steps; this allows us to control the decay of 
coefficients and, in such a way, to prove that v) < oo- 

In the next subsection, we handle the general case, in particular, 

• We give the generic calculation, which allows us to estimate the convolution 

{9t-sags){x)-, 

• We estimate the convolutions y), k >2; 

• We change the estimation procedure after ko steps, where ko is properly chosen, 

and estimate y), i > 1. 

The change of the estimation procedure could be unnecessary if we would know 
thatpt(x, y) possesses a more regular upper bound than (2.12). In this case, we obtain 
a sufficient control on the coefficients Ck, k > 1. This is exactly the case under the 
conditions of Proposition 4. 


3.3 Representation lemma, generic calculation, and estimation of convolutions 

Lemma 2. The function pf{x, y) given by (2.6) satisfies Eq. (3.1). 

Proof. In the case when is a symmetric stable-like process and w € Sk with 
respect to the transition probability density of X, the sketch of the proof is given 
in [22]. In the general case, the proof is the same; in order to make the presentation 
self-contained, we present it below. Using the equality 

e^* = [ e^*~^-^dAs + 1, 

the strong Markov property of X, and the additivity of At (cf. (2.2)), we write 


T,^f{x) = E^[f{Xt)e^^] 
= E^f{Xt) + 

= E^f{Xt)+E^ 

= E^f{Xt) + E"^ 



/‘ e ^» [fiXt-s)e 
Jo 

f Ttsf{Xs)dAs. 
Jo 



Observe that for / € we have 


f{Xs)dAs= f f f{y)ps{x,y)w{dy)ds. (3.11) 

Jo Jw^ 

Indeed, since xt = xt ~ Xf with xt given by (2.5) is the characteristic of At, 
Eq. (3.11) holds for a finite linear combination of indicators. Approximating / S 
by such linear combinations and passing to the limit, we get (3. 1 1). □ 
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For 9 e [0,1], put 

gt,d{x) := gt{9x), (3.12) 

where gt{x) is defined in (3.7), and 

:= p^~^^h{iyps), V > 0, (3.13) 

where h is the volume function (cf. (2.14)) appearing in condition (2.15). Lemma 
below gives the generic calculation, needed for the proof of Theorem 3. 

Lemma 3. For 9 £ (0,1), we have 

{gt-sags){x) < - s) + (j)(^i_g'jb{s)]gt,e{x), x e , 0 < s < t < 1, 

(3.14) 

where C > Q is some constant, independent of 9, and b > 0 comes from the definition 
ofgt, see (3.7). 

Proof. Take 9 £ (0,1). Since by definition the function pt is decreasing, we have 

\\x - z\\pt-s + \\z\\ps > ||a;||pt, 

which implies 

{gt-sags){x) 

< pt_,p^ [ [fup{{z - x)pt-s)fup{{y - z)ps)f~^'’\w\{dz). 

By integration by parts we derive, using that pt is monotone decreasing, that 

f [/up({a: - z)pt-s)fnp{{z - y)ps)] ^~^\w\{dz) 

< Pt /2 [ P^/up ® {{z - y)Ps)\rzj\{dz) 

pCO 

< Cip^p^ / \vu\{z : > e-'’}e-'^dv 

Jo 

poo 

= (1 - 9)bcip]^p^ / \w\{z : ||z - y|| < v/ps}e-^'^^-'^'>^dv 

JO 

pOO 

< (1 - 9)bcip'^p: / 

= c,p7p:+^h{b(i - 9)ps) 

= cip7A{i-9){s)- (3.15) 

Similar estimate holds true for s > |, which finishes the proof of (3.14). □ 

Take a sequence {9k)k>i such that 

Pi = l, Pfc>0, 9k-i>9k, k>2. (3.16) 


Let 


kn := 


n 

aC, 


(3.17) 
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where ( is the parameter appearing in (2.15). Define 

K := min{6(0j_i — Oj), 1 < J < fco}, 

F{t) := f (l)^{r)dr, 

Jo 

and 
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(3.18) 

(3.19) 


:= 


\gt,dkix)F'" (t), l<k<ko, 

I fc > fco, 


(3.20) 


where gt,e{x) is defined in (3.12). 

Finally, define inductively the sequence of measures 

'R^^\dw) -.= Qt{dw) if fc = 1, 


:= (2F(f)) ^ f f - s) + (j>K{s)]Qt-s{dw - ^^du)ds 

Jo JR" 

(3.21) 


if fc > 2. Since {Qt)t>o is the family of sub-probability measures (see Theorem 2), 
we have 


< (2F(f)) ' [ [(/)«(f-s) + <^Js)]Qt_«(M")Q,(R")ds<l, 

and we can see by induction that (R”) < 1, f G [0,1], for all fc > 2. 

Lemma 4. For fc > 2 we have 

\p^\x,y)\<Ck{g't'> a;,j/GK", fe (0,1], (3.22) 

(k)\ (k) 

where the sequence {gl ^)k>i is given by (3.20), ^ is defined in (3.21), k > 2, and 

for fc > fco, the constants can be expressed as 

Ck = 


where M, C > 0 are some constants. 

Proof. We use induction. Rewrite the upper estimate on pt(a;, y) in the form (3.5). 
For k = 2we get, using (3.5) and (3.15), the following estimates: 

rt p r p 1 


Pi 




y)\<Cf / / g^'l{z - X - wi)g[^y{y - z - W 2 )\Midz) 

Jo JR2" .JR" 


• Qt-s{dwi)Q s{dw2)ds 
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< C 2 



9t,e2{x 



- s) + (/)«(s)] 



{dw 


u)Qs{du)ds 


= 2C2{gf 


(3.23) 


where Ci > 0 comes from (3.5), and in the third line from below we used that by the 
definition of k and monotonicity of (j>i, in v. 


t e (0,1]. 

Suppose that (3.22) holds for some 2 < fc < fcg. Then 


\pf(-'^+^\x,y)\ <2^ ^CkCi f f {f]}^*Qt-s)iz - x) 

Jo Jr^ 

■ ( 5 ?^ *‘Jii'"^){y - z)dzds 

= 2^-^CkCi [ [ [ {y-x-wi- W2) 

Jo Jr'^ Jr'^ 

■ Qt-s{dwi)'36^\dw2)ds. (3.24) 

By the same argument as those used in the proof of Lemma 3, we have 

< {gt-sfi^o 9s,9k){.x)F’^~^ {t) 

< Ck+l9t,9k+Ax)F'"~^ {t) [(t)b(9t,_^-9k){t - s) + ct)b(g^_^_g^j{s)] 

= Cfe+i(F(f)) s) + (/i«(s)]pp+^^(a;). 

Substituting this estimate into (3.27), performing the change of variables and normal¬ 
izing, we get (3.22) for 2 < k < ko- 

Take cq > 0. Note that for some ci > 0, we have coPt < Pcit, t € (0,1]. Then, 
by (2.15), 

f Pt~^^h{coPt)dt <C2 [ PcfHPcit)dt <C3 f Pt^^h{pt)dt < . 

Jo Jo Jo 

Therefore, taking ko as in (3.17), we get 

p"F'=«(f) < < C6, te [0,1]. (3.25) 

In such a way, on the {ko + l)-th step, we obtain 


{gt°sngi^^){x) <ce-^^’^oPdM f |tn|(dz) 

Jr-^ 

pOO 

= / |zz7|{z: psb{l — 6ko)\\z — x\\ < r}e~'^dr 

Jo 
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< cg[^°^^\x)(j)^{s)F-^{t) 

(cf. (3.15)), where in the last line we used the inequality k < 6(1 — Oka) and the 
monotonicity of (j)^ in ly. Using this estimate, we derive 

<CkoCi f f f {g^'l°]ng^^'>){y-x-Wi-W2) 

Jo JR^ JR-^ 

■ QsidwO'Jil’^g {dw2)ds 

< 2cC,Ck, • “+'’)(!/ - x). (3.26) 

Then (3.22) follows by induction. Indeed, assume that (3.22) holds for k = ko + i—1. 
For f > 2 we get 

ogi^'’){x) < cF^~'^{t)e~'^'^'’oPdMfj)^{s) = cF-'^{t)g[^°~^^\x)(j)^{s). 

Therefore, 


I 0(^0+-^) 

\Pt 


{x,y)\ < (2Cic)^-iCiCfe„ r [ 

Jo Jr^ 

■ {gi^^ *Qs){y - z)dzds 
= Cfe„(2Cic)'(gf“+') * - x). 


){z-x) 


(3.27) 

□ 


Remark 4. As we observed in the proof, the estimation procedure depends on con¬ 
dition HI, which guarantees the existence of the number fcp such that (3.25) holds. 
In general, without HI we cannot guarantee the existence of such a number, which is 
crucial in our approach. For example, suppose that ps x s“^ for small s, and take the 
measure w such that 

7 - 5 -’ »'e(o,i]. 

In r 

By the Tauberian theorem, we have 6(A) x [Aln^A]“^ for large A. Therefore, 
4>v{t) ^ |lnf|“^ as f —>■ 0, and thus the integral F{t) diverges. Nevertheless, as¬ 
sumption HI can be dropped, if the function (a;, y) possesses a more precise upper 
bound. We discuss this question later in Section 4. 


3.4 Proof of statement c) 

From (3.27) we get for all x,y € K", 

\pf^'^°+^\x,y)\<M{CF{t)Y, e>l, (3.28) 

where M = Ckg and C = 2CiC. Without loss of generality, assume that C > 1. 
Since F{t) —>■ 0 as f —)• 0, there exists to > 0, such that 

CF(f)<l/2, fe(0,fo]- (3.29) 

Thus, for t € (0, to], the series (3.2) converges absolutely and is the solution to (3.1). 
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Let us show that the integral equation (3.1) possesses a unique solution in the 
class of functions {f{t, x,y) > 0, t € (0, to], x,y G K"}, such that 


/ f{t,x,y)dy<c, t G {0,to], x G . 
JR" 


(3.30) 


Then the series (3.2) is a unique representation of the Feynman-Kac kernel pf^ix, y) 
fortGiO,to],x,yGR-. 

Suppose that there are two solutions pf'^{x,y) and p'f'^{x,y) to (3.1). Put 
Pt{x,y) ■■= \pf'^{x,y)-pf'^{x,y)\ and r;t(x) := pf-{x,y)dy. Then, hy 
(3.1) we have 


By induction we get 


vt{x) < / / pt-six,z)vs{z)zu{dz)ds. (3.31) 

Jo JR" 

Vtix) < f f pf^^{x,z)vsiz)w(dz)ds. (3.32) 

Jo JR" 

Note that there exists c > 0 such that (x^ y) c for all t G (0, fp], x,y G M” 

(cf. (3.26)). In such a way, by the hniteness of measure w, we get 

Vt{x) < Cl / Vs{z)w{dz)ds <C 2 Vsds, (3.33) 

Jo JR" Jo 

whereug := supj,gjjn t;s(z). Taking sup^g^n in the left-hand side of (3.33), we derive 

vt<C2 [ Vsds, t G {Q,to]. (3.34) 

Jo 

Applying the Gronwall-Bellman lemma, we derive vt = 0 for all t G (0,fo]- Thus, 
the solution to (3.1) is unique in the class of functions 

{fit,x,y) >0,tG (0,fo], x,y G K”} 


satisfying (3.30). 

Estimating series (3.2) from above, we get an upper bound in (2.16) with /up of 
the form (2.10) and 

Jlt(dw) = cq y^ c^Jlt^^(dw), 

k>i 

with some c G (0,1) and the normalizing constant co > 0 chosen so that !Kt(R") < 1 
for all t G (0, fo]- 

For the lower bound, observe that by (3.20) we have 

|pf(x,t/)| <C(fco)pr^(f), 2<k<ko. (3.35) 

By (3.28) and (3.29) we get 

^ pf ix,y)< 2MCF{t), f G (0, fo], 

t>i 
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which, together with (3.35) and the observation that pt is decreasing, yields the esti¬ 
mate 





(3.36) 


where Co > 0 is some constant. Therefore, choosing to small enough, we have by 
the lower bound in (2.12) the inequalities 


ptix,y) > Ptfiow{{y - x)pt) - CoF{t)p^ 

>cptf\ow{{y-x)pt), te{0,to]. 


(3.37) 


□ 


4 Proof of Proposition 4 

Since the proof of the proposition follows with minor changes from the proof of 
the upper estimate in [22, Thm. 3.3], we only sketch the argument. For {t,x,y) € 
(0,fo] X M" X R”, put 



Io{t,x,y) := gt{x - y), Ik{t,x,y) = 


gt-s{x - z)Ik-i{s,z,y)w{dz)ds. 


By the same argument as in [22], we can get 


Ikit,x,y)\ < c'"Qt{y-x), k>l,t£ (0,fo], 


where c G (0,1) is some constant. Thus, for /c > 1, we have 

\p?’"{x,y)\ <c^3tix-y), a;,y G R”, f G (0,fo]- (4.1) 


This proves the convergence of the series (3.2) and the upper estimate (2.20). □ 


Remark 5. Let us briefly discuss the crucial difference between the proofs of Theo¬ 


rem 3 and Proposition 4. We changed the procedure of estimation of pf^{x, y) after 


a certain step, which was possible due to (2.15). In the case when we have a single¬ 
kernel estimate for pt{x,y), for example, (2.19), we can drop condition (2.15). In 
fact, it is enough to require that w G Sk with respect to Qtiy — x). This happens 
because in the case of the single-kernel estimate of type (2.19), it is possible to show 


that the convolutions pf^{x, y) satisfy the upper bound (4.1) with c G (0,1), which 


implies the convergence of the series (3.2). 

5 Examples 

As one might observe, the scope of applicability of Theorem 3 heavily relies on the 
properties of the initial process X. To assure the existence of such a process, we 
applied Theorem 2. Below we give the examplesin which condition A1 is satisfied. 
Since conditions A2-A4 are easy to check, we may assume that the functions a{x) 
and m(x, u) are appropriate. We confine ourselves to the case when the measure p in 
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the generator of X is “discretized a-stable; up to the author’s knowledge, in this case 
the cotTesponding Feynman-Kac semigroup was not investigated. Examples below 
illustrate that our approach is applicable also in the situation when the “Levy-type 
measure” m{x, u)iJ.{du) related to the initial process X is not absolutely continuous 
with respect to the Lebesgue measure. 

Example 1. a) Consider a “discretized version” of an a-stable Levy measure in R". 
Let mk,v{dy) be the uniform distribution on a sphere Sfc,^ centered at 0 with radius 
2 -kv^ t) > 0, A: G Z. Consider the Levy measure 

OO 

d-{dy)= Yl (5.1) 

k— — oo 

where 0 < 7 < 2v. In [17], it is shown that for such a Levy measure condition A1 is 
satished, and 

fe(0,l], (5.2) 

where a = 'y/v. 

Take some functions a(-) : K" —>■ K and a non-negative bounded function m(-, •) 
dehned on K" x K" satisfying assumptions A2-A4. By Theorem 2 the operator of 
the form (2.7) with /r, a{x), and m{x, u) as before can be extended to the generator 
of a Feller process X that admits the transition density pt{x, y) satisfying (2.12). 

Let 07 be a hnite Borel measure, and let h be its volume function, see (2.14). Let 
us show that if the inequality 

[ ^ ^ (0, 1], (5.3) 

Jo ^ 

for some C > 0^ then we have (2.15). Using (5.2), changing variables, and applying 
the Fubini theorem, we derive 



Denote by I{t) the right-hand side in this expression. Applying (5.3), we get 

poo 

I{t) < Cl / 

Jo 

In particular, if h{v) < cv^, d > n — a, then (5.3) holds. 

Thus, by Theorem 3, the Feynman-Kac semigroup {T^)t>o is well dehned, and 
the kernel pf{x,y) satishes (2.16) with some constants oi, 1 < i < 4, and some 
family of sub-probability measures 

b) Consider now the one-dimensional situation. In this case, the Levy measure p 
from (5.1) is just 

OO 

p{dy)= Y 2 "^(V..(dt/)+ 5 _ 2 -n.(dt/)). 


n— — oo 


(5.4) 
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Let X be a Levy process with characteristic exponent 


V '(0 := 


/ (l — C0S(^M 


))n{du). 


In [20] we show that if 1 < a < 2, then the transition probability density pt{x, y) of 
X, Xq = X, is continuous in (f, x, y) G (0, oo) x K. x K. and admits the following 
upper bound; 


Pt{x,y)<ct ^/°{l + \y-x\/f^°‘) f G (0,1], x,t/G K. (5.5) 


Note that the right-hand side of (5.5) is of the form (2.17) with d = 0. Thus, the 
conditions of Proposition 4 are satished, and we can construct the Feynman-Kac 
semigroup for the related functional At and the transition density pt{x,y), and get 
the upper bound for the function pf{x, y) with pt x t G (0,1]. 

To end this example, we remark that it is still possible to construct the upper 
bound for such pt{x, y) for a G (0,1) of the form but the function 

/ in this upper bound might not be integrable; see [20] for details. Note that the upper 
bound (5.5) is non-integrable in R” for n >2. 

Example 2. Consider the Levy measure 



'OO 


tAirv)r ^ drpo{dv), aG(0,2), (5.6) 


where a G (0, 2), po is a hnite symmetric non-degenerate (that is, not concentrated 
on a linear subspace of R") measure on the unit sphere S" in R". Suppose that there 
exists d > 0 such that for small r we have 


1^0 {B{x, r)) < Cr‘^, ||a:|| = 1. 


For d + a > n, it is shown in [6] that the corresponding Levy process X, Xq = x, 
admits the transition probability density pt{x, y), which satishes 


Pt{x,y)<ct ||y-a;||f f>0, x,yGR". (5.7) 


In the forthcoming paper [7], we construct a class of Levy-type processes that ad¬ 
mit the transition densities bounded from above by the left-hand side of (5.7). Thus, 
taking vu G we may apply Proposition 4. 
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A Appendix 


Proof of Lemma 1. We follow the idea of the proof of [4, Lemma 11]. Without loss 
of generality, assume that zu is non-negative. Suppose first that w £ 3C„_q. Using 
integration by parts and the Fubini theorem, we get 

n 0 s{a: - y)w{dy)ds 

Jo 

= f s~'^^°'zu{y : ||tc — y|| < 

Jo 


+ / s 


—njo 


1/a \ '^+' 


= a l-f 


\\x—y\\>s''-/°‘ V lU 


7{dy)ds 


a 


d + 2a — n 


w{y : llx - y\\ < u} 


,,n+l—a 


aid + a) d+ 2 a 

-t 


w{y : ||a;-y|| < v} 


d + 2a — n 




v‘ 


d+l+a 


dv 


dv. 


(A.l) 


Since va £ Xn,a, the first term tends to 0 as t ^ 0. Further, since d > n — a and the 
measure va is finite, we have 


d+2a-Ti 

t “ sup 


^{y ■ Ik - y|| < ■y} 


V 


d+l+a 


dv —^ 0 , t —^ 0 . 


Let us show that 


d+ 2 ^ p w{y ■. ||x-y||<t;} 
supf “ / -- dv. 


nd-t-l+Q: 


/tl/- 

Let Kq = Ko{t) := + 1; note that 1 as f —>■ 0. We have 

rf_e+ 2 ^-n P w{y : Ik-y||<'(;} 


(A.2) 


/(l/o 


V 


d+l+a 


-dv 




Kq /\ (d—n+2a)/a r ii II ^ 1 

sri^y I : Ik - y|l < 


kj 




nn+l-a 


dv. 


Since d>n — a,we have ^ n+ 2 a)/a < oo. Since zu £ X„^a, we have 

'•(^ + Vi ' . Ila- _ y|| < .y} 


max sup 

\<k<KQ{t) X Jkt^/^ 


,,n+l—a 


-dv —5- 0 , f —)• 0 . 


Thus, we arrive at (A.2). This proves that (2.18) implies that zz £ Sk with respect to 

9tiy-x). 

The converse is straightforward. □ 
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